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ABSTRACT 

Context. Most diffuse clouds are only known as kinematic features in absorption spectra, but those with appreciable Ht content may 
be visible in the emission of such small molecules as CH, OH, and CO. 

Aims. We interpret in greater detail the extensive observations of l2 CO emission from diffuse gas seen around the archetypical line of 
sight to £ Oph. 

Methods. The 12 CO emission is imaged in position and position-velocity space, analyzed statistically, and then compared with maps 
of total reddening Eg_ v and with models of the C + - CO transition in H 2 -bearing diffuse clouds. 

Results. Around f Oph, 12 CO emission appears in two distinct intervals of reddening centered near Eg_ v ss 0.4 and 0.65 mag, of 
which Si 0.2 mag is background material. Within either interval, the integrated 12 CO intensity varies up to 6-12 K km s _1 , compared to 
1 .5 K km s toward f Oph. Nearly 80% of the individual profiles have velocity dispersions cr v < 0.6 km s _I , which are subsonic at the 
kinetic temperature derived from H 2 toward ( Oph, 55 K. Partly as a result, I2 CO emission exposes the internal, turbulent, supersonic 
(1-3 kms~') gas flows with especial clarity in the cores of strong lines. The flows are manifested as resolved velocity gradients in 
narrow, subsonically-broadened line cores. 

Conclusions. The scatter between N(CO) and Eb-v in global, CO absorption line surveys toward bright stars is present in the gas seen 
around £ Oph, reflecting the extreme sensitivity of N( 12 CO) to ambient conditions. The two-component nature of the optical absorption 
toward f Oph is coincidental and the star is occulted by a single body of gas with a complex internal structure, not by two distinct 
clouds. The very bright 12 CO lines in diffuse gas arise at N(H 2 ) ~ 10 21 cirr 2 in regions of modest density n(H) as 200 - 500 cirr 3 
and somewhat more complete C + -CO conversion. Given the variety of structure in the foreground gas, it is apparent that only large 
surveys of absorption sightlines can hope to capture the intrinsic behavior of diffuse gas. 
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1. Introduction 

The line of sight to the nearby (140-160 pc) runaway 09.5 V star 
HD 149757, f Oph, has served as the archetype for detailed ob- 
servational studies of the internal composition of diffuse (Ay^ 1 
mag) clouds (Herbig 1968, Mort on) [1975[|, for opt ical/uv de- 
tection of new molecules in diffuse gas (Maier et al. 200 1\ and 



for theoretical models of molecular gas in diffuse clouds ( Black 



& Dalgarno 


[1977| 


et al."ni996). 


The H 2 



are dense enough to host appreciable column densities of carbon 



monoxide, N( 12 CO) « 2.4 x 10 15 cm' 2 (Morton 1975 



|et~aT) [1982) |Lambert etal) [T994) [SonnentruckereTa 



Wannier 
2007) 1 



and these are readily detectable in mm-wave emission toward 
the star ( |Knapp & Jura] 1 1976) |Liszt)| 1979) |Langer et aT)|1987)) 



They were very partially mapped in CO emission ( Kopp et al. 
~ Ltszt| [1997) 1, as well as CH and OH ( )Cratcher|[l979)|Liszt 



1996 



1997). 



CO J = l-0 emission around £ Oph was imaged in much more 



complete fashion by ITachihara et al. (2000), who focused their 



discussion on the properties of the nearby dark cloud complex 
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L204 seen several degrees to the galactic South. L204 is clearly 
outlined against th e Ha emission from the ionized gas in the 



star's H II region (|Gaustad et al.||200l|). In this work, the 12 CO 



datacube from Tachihar a et al.| ( |2000) > is employed to study the 
diffuse gas at Ay « 1 mag seen nearer the star. We scrutinize the 
entire CO image of the absorption-line host whose overall prop- 
erties have so often been inferred from one microscopic absorp- 
tion sightline toward the star, and we inquire to what extent that 
line of sight faithfully represents the host gas. Moreover, large- 
scale maps of reddening and extinction have become available at 
comparable resolution (though only along the entire line of sight, 
see |Schlegel et al.| ( fl998] ) and Dob ashi et aL) ( |2005| )), and we em- 
ploy these to control against possible confusion between diffuse 
and darker sightlines, a source of concern given the strong CO 
lines we see. 

The plan of this work is as follows. Section 2 summarizes 
what is known observationally of the line of sight toward the star 
and describes the pre-existing H I, CO, Ay and Eg v datasets 
which are discussed here. Section 3 discusses the appearance of 
the sky around ( Oph in terms of the statistics of 12 CO emission 
and reddening. Section 4 discusses CO profiles and linewidths in 
terms of the turbulent flows which are prominent in the emission 
profiles. Section 5 discusses physical conditions in the CO and 
H2-bearing host gas, especially the regions of extremely bright 



Liszt, Pety and Tachihara: CO, reddening and turbulence around ( Oph 



23 




L — 6.0<v<9,0 km s 1 Iflflj 

a=1 6h54m 

GALACTIC LONGITUDE 



Fig. 1. Limiting reddening (Schlegel et al. 



1998) and 12 CO J=l-0 emission (Tachihara et al. 



2000) in the vicinity of f Oph. The 



12 CO emission contours are shown at levels 1, 2, 4, 6, 8, 10, 12, 16, 20, 24 K km s" 1 and have been calculated over different velocity 
ranges in the upper right and lower left portions of the map, as indicated at the map corners. The dark cloud complex to the South 
of £ Oph is known as L204 and the gas at the upper right, around f Oph, as L121. 



(1 1-12 K) CO emission. Sect. 6 is a summary and Sect. 7 (avail- 
able online) discusses the relationship between the extinction 
and reddening measurements over the £ Oph field and presents 
some additional views of the 12 CO observations. 



2. Observations 



2.1. Carbon monoxide 



The datacube of Tachihara et al. (2000) comprises nearly 1 1,000 
spectra from the NANTEN telescope with a beamwidth HPBW 
= 2.7' on a 4' grid in galactic coordinates. The spectra have 0.1 
kms" 1 resolution and the single-channel rms at this resolution, 
0.5 K, is relatively high compared to that in the small numbers 
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Fig. 2. Integrated H I line brightness from the LDSS survey 
of Hartm ann & Burton] (JT997J at 35' resolution, converted to 
equivalent reddening E B -v= N(H I)/5.8 x 10 21 cm 2 , N(H I) 
=1.823 x 10 18 cm- 2 j Tb(H I dv)/(K km s~'). 



of demonstration spectra typically shown in earlier work (Liszt 
1997) . Statisically significant detections of the CO require Wco 
^ 1 kms' 1 , where Wco is the integrated intensity 

Toward the star, we show the profile of |Liszt| ( |1997| > from 
the then-NRAO Kitt Peak 12m telescope at 1' (HPBW) spatial 
resolution and 0.12 kms~' spectral resolution (see Sect. 3) . To 
ensure that this profile is compatible with those from NANTEN 
we recently used the ARO 12m Kitt Peak telescope to re-observe 
several positions having comparatively strong emission in the 
NANTEN datacube. The 12m spectra agree with the NANTEN 
data to better than 5%, a remarkable coincidence considering the 
difference in hardware and spatial resolution. 

2.2. H I 

To help distinguish between foreground and background mate- 
rial, or atomic and molecular gases, we employed the H I pro- 



files from the Leiden-Dwingeloo all-sky H I survey (Hartmann 
& Burton 1997[ ). These data have 35' resolution on a 0.5° grid 
in galactic coordinates. 



2.3. Reddening and extinction 

The reddening toward £ Oph is known to be Eb v = 0.32 mag 
( |Morton||l975| >, but maps of the foreground extinction over only 
the first 140 pc are unavailable. Instead, we employ the limiting 
reddening (from here to infinity) from the work of Schlegel et al. 

v , with a spatial resolution of 6.1', pub- 



1998] l, denoted by E™_ 

lished on a 2.5' grid; the stated global rms error of this dataset 
is 16% (a percentage at each pixel). Their values are based on 
a determination of the dust column density estimated from the 
IRAS 100 micron flux adopting the temperature variation de- 
rived from COBE/DIRBE 240 micron data. The minimum lim- 
iting reddening in the region, approximately 0.23 mag (Fig. 1), 



is likely hosted in atomic gas over long paths, as opposed to the 
more-localized diffuse and dark clouds of interest here. 

A comparison between the limiting reddening of Schlegel 
et al. ( 1998 1 and maps of the limiting extinction A^ computed 



from star counts by Dobashi et al. ( 2005 1 can be found in Sect. 
7 (online only). 



2.4. Some general conditions in the gas along the line of 
sight to £ Oph 



In front of the star, E B -v = 0-32 mag, N(H I) = 5.2 x 10 



20 



N(H 2 ) = 4.5 x 10 20 cm" 2 (Savage e tai]|l977ll, N( 12 CO) 

. „ic i L _ 1 I . i . — ' ' . I I . J ' .1 I _ 



10 15 cm' 2 {Morton 
|etal.||2007LN(Cl 



1975 

T 



Wannier et all 11982 



cm , 
2.4 x 



Sonnentrucker 



MrCQl jMortonnT975| land N(C + ) > 3 x 



10" cm ^Cardelli et"aH|l993| i7The" mean kinetic temperature 
of the molecular gas inferred from H 2 absorption is 54 K ( Savage 
etaLl[T977] l. 

Toward the star the limiting extinction from the work of 
Schlegel et al.|(|1998|l is Eg_ v = 0.55 mag, so that the background 



reddening is approximately 0.23 mag. 

The distance to the occulting material is generally taken to 
be very close to that of the star, just outside the nearer edge of 
the star's H II region ( [Wood et aL] |2005| >. At a distance of 140 
pc, 1' corresponds to 0.041 pc and 1° to 2.44 pc. 



3. The sky around C Oph viewed in reddening and 
CO emission 

Figure 1 is a composite image of the limiting reddening Eg v 
and integrated 12 CO J=l-0 emission, an updated version of Fig. 
1 



Tachihara et al. ( 2000|l. The gray-scale underlay is the 
limiting reddening ( Schlegel et al.| 1998] > normalized to white 



at the minimum value seen over the region, Eg v = 0.23 mag. 
Superposed on the reddening map, the integrated I2 CO intensity 
Wco has been calculated separately over the diffuse/translucent 
northwest region, referred to as L121, and the translucent/dark 
southeast region, L204. For the diffuse gas of L121 at upper 
right, the red and blue contours correspond to the velocity ranges 
above and below v = -0.25 km s' 1 as indicated in the upper right 
corner: this division corresponds to the natural separation be- 
tween the two components of CO, CH and OH emission found 
toward and around £ Oph, as shown in the profiles and position- 
velocity maps of |Liszt| ([T997). For the darker gas of L204 to the 
southeast Wco was calculated over three intervals but nearly all 
of the emission from L204 occurs in the velocity interval 2-6 
kms -1 represented by the white contours (see Fig. 5). The noise 
level is somewhat larger for these contours, because Wco was 
calculated over a somewhat broader interval. 

The strong CO emission associated with L204 often follows 
the ridge lines of the extinction with something of a setback (red- 
dening without apparent CO) in the direction of the star. This 
is consistent with the edge-on geometry for L204 described by 
Tachi hara et al.| (|2000). Nearly all of the CO emission in the 
L204 dark cloud complex is found in the interval V/ sr = 2-6 
kms" 1 and so does not overlap that of the more diffuse gas to 
the north, seen at v/ sf % 2 kms 1 . The extent to which the two 
regions are at rest in the directions joining them, and might par- 
take of the optical pumping excitation mechanism for diffuse gas 
described by |Wannier et al.| ( |1997| l is unknown. |Tachihara et al.| 
(2000) assumed that the two clouds were co-moving at the edge 
of the H II region around ( Oph in order to discuss the energetics 
of the gas. 
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Fig. 3. Integrated intensity of the CO line Wco vs. limiting red- 
dening y . Top: for the diffuse region LI 21 at upper right in 
Fig. 1; bottom: for the darker L204 gas at lower left in Fig. 1. 
In each case the data are shown twice, at right as individual data 
points, at left binned and contoured at logarithmic (factor two) 
intervals: the contour representation of the L121 data is super- 
posed at bottom. Also, in each case the data are divided into 
regions of lower and higher extinction with the separation oc- 
curing at the mean Eg v over datapoints which lack statistically 
significant detections of CO emission (0.46 mag and 0.67 mag at 
top and bottom, respectively). The locus of the sightline toward 
( Oph in the B-region of the upper panel is indicated by a (red) 
cross. Mean 12 CO emission profiles are shown in Fig. 5. 



To the North of L204, on either side of the star is the ( Oph 



diffuse cloud, described as the L121 complex by Tachihara et al. 
( 2000| l. Most of the emission from that gas occurs at -2.4 km s -1 
< V/ sr < 2.2 km s . The red-shifted kinematic component at -0.2 
kms 1 < v/ SJ - < 2.2 km s _1 which is seen in CO emission and in 
many species in absorption toward the star, is very much con- 
fined to the northern and northwestern edges of the broader dis- 
tribution of blue-shifted gas. Its separate identity as a second 
cloud is somewhat marginal in the emission maps and hardly 
supported by a more detailed examination of the turbulent gas 
kinematics (also see Fig. 7). 

The diffuse gas in the L121 complex around ( Oph is sep- 
arated from the L204 dark cloud by an extended trough in 
the reddening whose overall mean value is Eg v « 0.43 mag. 
Immediately below the star is a pronounced minimum whose 
mean is Eg v * 0.34 mag and whose absolute minimum is Eg v 
= 0.29 mag. The overall impression is of a cylindrical shell ge- 
ometry and perhaps a separate, more circular shell of radius * 
1° around the star. However, a yet-larger map of the extinction 
shows that L204 is part of a much larger ridge. 



3.1. Separating foreground and background reddening 

As noted above, we wish to understand the relationships between 
Wco and reddening, but we only have maps of the reddening 



which include all of the unrelated background and foreground 
gas. Immediately around £ Oph Eg v = 0.55 mag, compared to 
Eb-v = 0.32 mag in front of the star, so some 0.55-0.32 mag 
= 0.23 mag of reddening occurs behind the star. This also cor- 
responds well to the absolute minimum in Fig. 1, Eg v = 0.23 
mag found some 1.5° North of the star. 

To estimate the amount of extraneous foreground material, 
we note that this is expected to be in atomic form and, in absorp 



tion against the star, N(H I) = 5.25 X 10 2U ( |Savage etal.||1977 i 
corresponding to Eb-v 55 0.1 mag. This leaves Eb-v ^ 0.22 mag 
of foreground material associated with the pure H2 component 
and gives the general idea that perhaps as much as 0.32 mag 
should be subtracted from the map of Eg v to infer the local 
reddening intrinsic to the H2-bearing gas. However, this is prob- 
ably an overestimate of the required correction because some of 
the foreground atomic gas is in the vicinity of LI 21 and L204, 
providing shielding. 

Figure 2 is a map of the integrated H I brightness, 
scaled to the mean relationship between hydrogen column 
density and extinction, that is, the map shows 1.823 x 
10 18 cirT 2 /T B (HI)dv/5.8 x 10 21 cm^mag- 1 where the veloc- 
ity integral is in units of K km s . Several of the local reddening 
peaks which lack CO emission in Fig. 1 even at Eg v = 0.65 mag 
are present as peaks in Fig. 2, for instance at (l,b) = (7.5°,21.5°) 
and (4°,20°). Some of this gas must be indigenous to the region 
of interest, but the extinction associated with the H I peaks is 
only w 0.06 mag, judging from the peak levels in Fig. 2 (0.28 
mag) compared to the nearby background level (0.22 mag). 

The minima in Fig. 1 and Fig. 2 and the difference between 
the foreground and limiting reddening toward the star consis- 
tently imply a background reddening contribution of * 0.2 mag 
over the region of interest. 



3.2. Quantitative relationship between CO emission and 
reddening 

To quantify the relationship between reddening and CO emis- 
sion, and to compare and contrast the diffuse and dark sightlines 
we divided the extent of Fig. 1 along a Northeast-Southwest 
diagonal in the trough of reddening between L121 and L204, 
along a line running from (l,b) = (3°,19.333°) to (10°,25°). Fig. 
3 shows the integrated CO intensities Wco an d limiting redden- 
ing Eg v for all points in both regions: the profile integral was 
taken over the range -2.4 kms -1 < v;„ < 2.2 kms -1 for the dif- 
fuse L121 gas shown in the top panel, and 2 kms _1 < v/ sr < 8 
kms" 1 for the darker L204 region in the lower portion of the 
map of Fig. 1 and the lower panel in Fig. 3. The rms noise in 
integrated intensity is 0.5-0.6 Kkms 1 . Profile integrals above 1 
K km s _1 generally represent real detections. 

Once we realized that the CO emission in the diffuse region 
was bimodal, as illustrated in the top panels of Fig. 3, we further 
sub-divided the diffuse gas into A and B portions corresponding 
to the two branches of the emission distribution. The A and B 
portions were separated at Eg v = 0.455 mag, which is the mean 
Eg v for those sightlines along which Wco < 1 K kms 1 and 
which therefore lack statistically significant detections of CO 
emission. The A-branch pixels have strong CO emission at Eg v 
substantially below the mean of those sightlines lacking CO 
emission at all. For the diffuse gas in L121 there is actually 
a substantial spatial segregation of the A and B portion pixels, 
with unweighted mean < (Lb) > = (6.0 o ±0.7°,23.8 o +0.7°) and 
< (l,b) > = (4.0 o ±0.9°,22.3 o ±1.0°) for the A and B portions, re- 
spectively. These centroids are on opposite sides of ( Oph and 
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Fig. 4. Cumulative probability distributions of the integrated in- 
tensity Wco (left) and peak line temperature (right) over the dif- 
fuse, northeastern L121 (at top) and dark southeast L204 sight- 
lines shown in Fig. 1. Only those sightlines at which Wco ^ IK 
kms -1 have been considered. 



separated by more than lcr + 1 cr in each coordinate. The red- 
shifted gas appears mostly in the higher-extinction B-portion 
while the blue-shifted gas appears more nearly in both the A 
and B-portions. 

In either emission branch, Wco varies widely over a rela- 
tively narrow range in Eg v , mimicing the extreme sensitivities 
of N(CO) which are seen in global absorption line surveys of sin- 



3.3. Statistics of the line brightness 

Figure 4 presents the distribution of integrated and peak bright- 
ness. The differences between the A (lower Eg_ v ) and B por- 
tions of the diffuse gas are somewhat more pronounced in the 
distribution of the line profile integral at left and somewhat less 
so in the peak, so that the line widths differ more than the line 
heights. Although the median brightness is at most only 3 K or 3 
K kms^'in the diffuse gas and the tail of the distribution seems 
very poorly populated above, say, 8 K peak brightness: paths 
which traverse LI 21 even in the short dimension have a very 
substantial chance of containing at least one such bright line, as 
discussed in Sect. 5 (see Fig. 7). 

It is not the intention here to discuss the dark gas, but it 
should be noted that the difference in mean brightness between 
the diffuse and dark regions are modest and correspond approx- 
imately to the differences in Eg_ v , thereby preserving the pos- 
sibility of a common CO-H2 conversion factor; the same is also 
true for the A and B portions of L121, see Fig. 5. This occurs 
despite the fact that most of the free gas-phase carbon is in C + in 
L121 (99% toward £ Oph) and in CO in the dark gas of L204 
(where N(CO) ~ 3 x 10 17 cirT 2 ) implying a difference in CO col- 
umn density and Wco/N(CO) by a factor of order 50. In Section 
5 we discuss the very different proportionality Wco <* N(CO) 
which is observed within the diffuse regime alone (Liszt 2007a). 



gle sightlines over widely-separated regions (see Liszt (2007a l) 



The widths of the two branches are comparable to the stated rms 
noise of the Eg v data (16%) and they are separated by 0.2 mag, 
which is of the same order as the background contribution: it 
is possible that they would more nearly coincide in Eb v if the 
background contribution were strongly variable. However, the 
B-branch at higher Eg v is more heavily populated at Wco > 5 
K kms -1 , suggesting that it actually is somewhat more strongly 
shielded, fostering a higher CO abundance and brightness: in this 
regime, Wco x N(CO) (ibid and see Sect. 6) and the conversion 
of free carbon from C + to CO occurs over a very narrow interval 
in Eb-v and/or NXH2). 

The distribution of brightness in the darker L204 region is 
not similarly bimodal at all Eg v and there is less rationale for 
a simple division into sub-portions. However when this is done, 
at Eg v = 0.67 mag corresponding again to the mean over pixels 
lacking statistically significant CO emission, the resultant lower- 
Eg v portion corresponds (in Eg v and Wco) to the entirety of 
the diffuse region. Consequently the lower-Eg v part of L204 is 
labelled L204 AB and the other L204 C. The behavior of W co 
with Eg v in the dark gas is complex but clearly bimodal for 
Wco ^ 12 Kkms -1 ; the strongest emission is by no means lim- 
ited to the darkest regions. Statistics of the brightness distribu- 
tion over L121 and L204 are shown in Fig. 4 and discussed in 
the following sections. 



3.4. An incidental bound on N(H 2 ) over the diffuse gas 

Because the line of sight to ( Oph occurs at such a high value of 
Eg v relative to the rest of L121 we infer that N(H2) is never very 
much larger in L121 than toward the star. For instance, if we take 
Eg v = 0.65 mag characteristic of the B-portion and subtract a 
background contribution 0.23 mag equal to that toward the star, 
the remaining gas column with Eg v = 0.65 - 0.23 = 0.42 mag 
corresponds to N(H) w 2.4 X 10 21 H-nuclei cirT 2 and N(H 2 ) 
< 1.2 x 10 21 cirT 2 . Conversely, because the line of sight toward 
the star has such a high value Eg v = 0.55 mag relative to the 
rest of the region, emission from the A-branch at Eg v = 0.4 
mag, probably arises in regions whose reddening and N(H2) are 
actually below those seen toward the star. 



4. Line profiles, linewidths and turbulent flows in 
diffuse gas 

4.1. CO line profiles and profile widths 

For a pure-H2 gas, the Doppler temperature equivalent to a 
given linewidth is T dopp = 43.6 K FWHM 2 or T dopp = 242 K 
cr v 2 where <x v is the velocity disperson in kms -1 . Typical ki- 
netic temperatures in CO-bearing diffuse gas are above 30 K 
( [Sonnentrucker et al.| [20071 l B urg h et [2007] pszt] |2007a) > 
and the mean kinetic temperature seen in surveys of H 2 absorp- 
tion is 70 - 80 K dRachford et al.[ [2002]|Savage et al.[[T^T7] >. CO 



profiles with FWHM ^ 1 kms _1 are subsonic at diffuse cloud 
temperatures and the purely-thermal velocity dispersions of 
CO molecules, 0.1 kms 1 at T K = 50 K, do not contribute im- 
portantly to the observed linewidths. 

Unweighted mean profiles formed over the sub-portions of 
the diffuse L121 and dark L204 regions are shown in Fig. 5. The 
mean profiles have linewidths which are supersonic, FWHM of 
typically 2-3 kms -1 but the individual sightlines in L121 typi- 
cally have subsonic linewidths; the FWHM of the two kinematic 
components seen toward the star at 1' resolution (shown in Fig. 
5 as the dark dashed line) are 0.6 kms 1 and 1.1 kms 1 , equiva- 
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Fig. 5. Mean spectra over the diffuse and dark sub-portions de- 
fined in Fig. 3: the profiles in the diffuse gas (L121) are those at 
lower velocity. The integrated brightnesses are Wco = 2.9 and 
4.0 Kkms-' for L121 and 2.3 and 7.4 K km s" 1 for L204. Also 
shown (dashed black line) is a spectrum at 1' resolution toward 
(, Oph from the 12m telescope. 



lent to Doppler temperatures of only 16 K and 53 K in a pure H2 
gas, compared to a measured temperature in H2 of 54 K as noted 
in Sect. 2.4. 

The distribution of velocity dispersions found over the LI 21 
region is shown in Fig. 6. To produce this figure we used the 
following windowing technique to measure the widths of spec- 
tra which might contain more than one kinematic component: 
select the overall velocity interval of the diffuse gas; find the 
peak channel; select those contiguous channels around this peak 
with temperatures above a noise threshold of 0.25 K; calculate 
the brightness-weighted velocity dispersion over those channels; 
mask off that portion of the profile; repeat until no channel above 
1 K remains unmasked. The dispersions measured en masse in 
this way are subject to overestimating the width in cases of un- 
recognizd blending, but they agree to within 10% for a series of 
test profiles which were fit with gaussian components (for in- 
stance, those shown in Fig. 8). 

The general properties of the gas in L121 have often been 
inferred from profiles seen along the single, microscopic absorp- 
tion line of sight toward the star and a sensitive CO profile to- 
ward ( Oph at 1' spatial resolution is also included in Fig. 5. 
How representative is it of the larger-scale distributions of the 
host gas? The stronger CO component at v; sr = -0.7 kms -1 to- 
ward £ Oph is one of the narrowest lines known in a diffuse 
cloud, with FWHM = 0.60 kms -1 . Broad consideration of this 
question is given in the discussion of Fig. 7 which displays most 
of the profiles seen in the L121 region, immediately following. 

4.2. Velocity gradients and turbulent flows 

The connections between the profiles seen at individual pixels 
and the mean profiles shown in Fig. 5 are the flows and velocity 
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Fig. 6. Statistics of line profile velocity dispersions in the dif- 
fuse gas. Large outer panel: cumulative probability distribution 
histograms of the measured velocity dispersion <x v shown sep- 
arately for the A and B portions of the diffuse gas (see Sect. 4 
of the text). The one-dimensional thermal velocity dispersions 
corresponding to kinetic temperatures of 20, 40 and 80 K in a gas 
of pure H2 are shown at upper left. Inset at lower right: variation 
of cr v with Wco- 



gradients in the host medium, i.e. the character of the turbulence. 



Although very elaborate analyses of line profile centroids (Pety 
|& Falgarone| 2003[> and fluctuations in t he wings of optically 



thick profiles ( Falgarone & Phillips 1996 1 have been used to in 



fer the properties of turbulence in denser neutral media, turbulent 
flows in the nearby diffuse gas of L121 are immediately visible 
in the line cores and (usually) spatially resolved into the shifts of 
individual, subsonic profiles. 

This is illustrated for £ Oph in Fig. 7, where we show a com- 
prehensive series of longitude-velocity diagrams spaced at one- 
pixel (4') intervals in galactic latitude, as labelled in the individ- 
ual panels. The peak line brightnesses (K) seen in the individual 
panels of Fig. 7 are labelled on the bars showing the color scal- 
ing in each panel. Although lines with peak brightnesses above 
8 K seem relatively rare in Fig. 4, they are common enough that 
peak brightnesses 7.5 K and higher appear in 40%, 18 of the 
44 individual panels of Fig. 7. As shown in Fig. 3, the most in- 
tense lines in the diffuse gas are by no means limited to the more 
heavily-reddened sightlines. 

Figure 7 shows that a description of the gas in terms of two 
identifiable foreground clouds at the absorption line or CO emis- 
sion line velocities toward f Oph is not appropriate. The line of 
sight to the star could just as well have occured behind any of 
the profiles exhibited in Fig. 7, leading to a wide variety of pos- 
sible interpretations of the intervening medium. Discussing the 
map in Fig. 1, we saw that the red-shifted component generally 
appears as something of a fringe at the northern edges of more 
broadly-distributed, negative-velocity gas. Viewed in position- 
velocity space in Fig. 7, the red-shifted gas is often seen as a pro- 
nounced kinematic excursion or wing, for instance, at b = 22.8° 
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or 23.4°. Further to the North, at b > 23.2°, the kinematic pat- 
tern undulates across the positive-velocity portion. At b= 23.4° 
the resolved velocity gradient spans the entire range of velocities 
in the diffuse gas: viewed from a different perspective, this same 
region might have been seen as a single, broad line. 

Figure 8 is an expanded view of the panel in Fig. 7 at b= 
22.6°, where we also show several included line profiles, their 
gaussian decomposition and the resultant FWHM. Except at the 
center of the diagram where there is a partially spatially -resolved 
velocity gradient, the FWHM are small enough to be subsonic 
for an H2 gas at typical diffuse cloud temperatures above 30 
K. Thus the observations show in detail how supersonic profiles 
might arise from the coincidental superposition or addition (for 
instance through beam-smearing) of velocity-shifted quiescent 
regions. 

As in Fig. 8 there are many highly-localized, relatively broad 
lines in 4' pixels, often joined to abrupt but spatially-resolved ve- 
locity gradients and reversals in narrow-lined gas. This suggests 
that broader lines are themselves composed of unresolved ve- 
locity gradients and it seems possible that any profile in Fig. 7 
with a width substantially above sonic is an unresolved gradient. 
For instance, compare the velocity span at b=22.87°, which is 
spatially resolved, with that at the edge of the emission region at 
b=23.33°. Recalling the extremely narrow blue-shifted CO com- 
ponent toward £ Oph and the low implied Doppler temperature, 
we wonder whether and at what scale profiles having linewidths 
that are subsonic at the higher temperatures expected for diffuse 
gas, say 60 K, might show velocity or spatial substructure. The 
geometry of the turbulent flows producing the patterns in Fig. 7 
will be considered in a subsequent paper. 

5. Physical conditions in the diffuse gas 

There are some 1000 sightlines with Wco ^ 1 K kms~' in 
the diffuse, Northwest portion of the map in Fig. 1. At a dis- 
tance of 140 pc, the H 2 -mass associated with these sightlines, 
parametrized by their average H2 column density <N(H2)>, is 
M = 430 M < N(H 2 )/1 x 10 21 cirr 2 > where typical values of 
N(H 2 ) along the L121 sightlines are 0.5 - 1.2 x 10 21 cirr 2 (see 
Sect. 2.4 and 3.4). This may be compared with the value 520 Mq 
given by Tachihara et al. ( 2000 1 based on a CO-H2 conversion 
factor N(H 2 )/W C o = 1.56 xKF cm^/XKkms 1 ). 

Nonetheless, deriving the physical parameters of host dif- 
fuse gas from CO profiles is challenging. Strong fractionation 
of the carbon isotopes (Liszt & Lucas| 1 1 9981 ISonnentrucker 



let al. 



[20571 |Burgh et al.| |2007| |Liszt| |2007a| ) causes the 



N( lz CO)/N( l3 CO) to vary in the range 15 - 150, thus making it 
impossible to derive the line excitation temperatures and optical 
depths, or the kinetic temperatures and column densities, under 
the usual assumption (valid in dark gas) that the relative abun- 
dances of 12 CO and 13 CO only reflect the general interstellar 
carbon isotope ratio. However, it is straightforward to show that 
the general properties of the CO observations, even including the 
rather unexpectedly bright lines, fit easily into the framework of 
diffuse gas at typical temperatures 30 - 60 K and modest density. 
The real underlying mysteries are the working of the poorly- 
understood polyatomic chemistry which forms the CO and other 
species at such modest densities ( |Liszt et aT]|2006| l, and the ori- 
gin of the turbulent flows which may power the chemistry. 

Figure 9 reports some results from the models which were 
used to interpret CO absorption line data by Liszt (2007a). The 
underlying physics are: to heat a uniform-density gas sphere 
by the photoelectric effect on small grains, as in the work of 
Wolf ire et al.| ( 1995 2003}, to calculate the ionization balance 



including grain-assisted neutralization of atomic ions (including 
most importantly the protons, ibid), to allow equilibrium of the 
self-shielding of H 2 formed on grains and the formation of CO 
by thermal electron-recombination of a fixed quantity of HCO + 
X(HCO + ) = N(H CO + )/N(H 2 ) = 2 x 10" 9 (the actual secular evo- 
lution is traced by |Liszt| ( 2007b[ >); and, to calculate the rotational 
excitation of CO assuming microturbulent radiative transfer with 
a linewidth determined by the local sound speed. Note that the 
interstellar radiation field in the models has not been increased 
above the mean interstellar value to account for the presence of 
the star and that all of the aspects shown here for the model re- 



sults are present in the absorption line data summarized by Liszt 
fl2007a| i. 

Model results are shown for just two fairly high densities 
n(H) = 256 cirr 3 and 512 cirr 3 . To form the graphs in Fig. 9, re- 
sults were derived by integrating along sightlines ranging across 
the faces of model spheres from center to edge. The graphs sum- 
marize results gleaned from models whose central column den- 
sities N(H) varied widely, so that the same value of N(H 2 ) might 
occur at different positions across the faces of models with dif- 
fering n(H) and N(H), and therefore have different N(CO) and 
Wco- For N(H 2 ) = 1 x 10 21 cirr 2 the models have typical sizes 
of 1.2-2.5 pc but the CO abundance is concentrated into smaller 
central portions of the host bodies owing to the chemistry of CO 
and that of H 2 , and the CO emission is more concentrated still 
owing to geometry and radiative pumping. 

As noted above N(H 2 ) = 0.5 - 1 X 10 21 cirr 2 in the CO- 
emitting regions around the star. As shown in panel c at the lower 
left, this is precisely the regime where carbon is about to recom- 
bine fully to CO at the quoted densities: the increase of N(CO) 
with N(H 2 ) is very rapid. Both N(CO) and Wco vary rapidly 
and have large scatter when plotted against N(H 2 ). Substantial 
CO column densities can accumulate in gas which is still rela- 
tively warm, 30-50 K, giving rise to 12 K lines as observed in 
the brightest profiles. 

The most important consideration is as shown at the upper 
right in panel b, W C o <* N(CO) for W C o ^ 10 K kms -1 , explic- 
itly independent of density and implicitly independent of N(H 2 ) 
and other cloud properties. This is a very general consequence 
of very sub-thermal excitation, as first shown by Goldreich & 



Kwan ( 19741 and does not require low optical depth. As shown 



in Fig. 4, some 80-90% of the diffuse cloud spectra have Wco 
< 5-6 K kms 1 and virtually all are below 10-12 K kms 1 , 
just in the regime characteristic of sightlines studied in uv and 
mm-wa ve CO absorption work generally, see Fig. 6 of |Liszt| 
(2007a 1. Therefore, the l2 CO brightness map in Fig. 1 is a map 
of N{ lL CO) in the diffuse gas and the same would be true for 
N( I3 CO) and the brightness of the 13 CO line. This is the one 
unambiguous result of mapping CO emission in any diffuse gas. 

The extreme sensitivity of Wco to N(H 2 ) therefore arises be- 
cause Wco 00 N(CO), so that a plot of Wco vs. N(H 2 ) is equiv- 
alent to plotting N(CO) against N(H 2 ). The net result is that al- 
though commonly-used values of the CO-H 2 conversion factor 
apply to some gas (as indicated in panel a of Fig. 9) the CO- 
H 2 conversion factor varies widely in diffuse gas and the actual 
N(CO)/N(H 2 ) ratio is small but very uncertain. A map of CO 
emission from diffuse gas is an image of the chemistry, not the 
mass distribution. 

Last, note that the J=l-0 CO line brightness is insensitive to 
density at fixed N(CO), indicating that other tracers are required 
to measure the local density when mm-wave emission profiles 
are analyzed. The J=2-l/J=l-0 line brightness ratios at lower 
right in Fig. 9 are not very sensitve to density, which explains 
why line brightness ratios 0.7 - 0.75 are indeed so commonly 
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Fig. 7. Longitude-velocity diagrams traversing the f Oph diffuse cloud at many latitudes within the Northwest portion of the region 
shown in Fig. 1. The vertical separation between panels is 4' (1 pixel) and the galactic latitude within each panel is labelled. The 
color scale in each panel are separately normalized within a range (Kelvins) shown on the individual color bars. 



observed in diffuse and translucent gas (Falgarone et al. 1998 6. Summary 



Pety et al.| |2008| l. The J=3-2/J=l-0 brightness ratio is a better 



indicator of density in the CO lines, but care must be taken to 
match the spatial resolution of the two lines. 



The line of sight to the nearby (140-160 pc) runaway 09.5V star 
£ Oph has for many years been used as an archetype for study- 
ing the properties of diffuse clouds in optical and uv absorption. 
Because the material has an appreciable molecular content, the 
host diffuse clouds can actually be imaged on the sky in space 
and radial velocity. Because the gas is well extended and com- 
paratively close, it provides an unusual opportunity for study of 
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V LSR [KM S" 1 ] 

Fig. 8. The longitude-velocity diagram at b = 22.6° from Fig. 7, 
expanded and decomposed into its constituent profiles at three 
positions (shown inset). Overlaid on each inset spectrum is a 
one-component gaussian fit whose FWHM is shown at the upper 
left. 



ing and turbulence around £ Oph 9 




Fig. 9. Models of CO formation and excitation. Results are 
shown for models having n(H) = n(H I) + 2n(H2) = 256 cm -3 
and 512 cm -3 . Upper left: Variation of integrated 12 CO bright- 
ness Wco with NXH2); the shaded line labelled "2" is a CO- 
H2 conversion factor 2 x 10 20 H2 crrT 2 (K kms" 1 )" 1 . Upper 
right: variation of Wco with N( 12 CO). Bottom left: variation 
of N( 12 CO) with N(H2). Lower right: comparison of I2 CO line 
brightnesses in the lowest two rotational transitions. 



diffuse gas and its interaction with its surroundings, including 
the star. 

We began by comparing the 12 CO J=l-0 emission line dat- 
acube of lTachiharalFaTl l pO&Ol l (HPBW 2.7' observations on a 4' 
= 0. 1 pc pixel grid) with a map of the limiting reddening Eg v 
from the work of Schle gel et al.| ( | 1998[ > having similar resolu- 
tion 6.1' (Fig. 1). The reddening in the L121 complex near and 
around f Oph ranges from Eg_ v = 0.23-0.75 mag and the inte- 
grated CO emission up to Wco = 12K kms , with 12 K peak 
temperatures, which are very bright lines indeed. Comparison 
of reddening of the star (0.32 mag) and through the Milky Way 
(Eg v = 0.55 mag) in the same direction, and comparison of H 
I seen toward the star in uv absorption and around the star in 
21cm emission (Fig. 2), suggests that as 0.2 mag of Eg v should 
be ascribed to unrelated background material. 

CO emission from diffuse L121 gas seen around £ Oph is 
bimodal in Eg v , clustering around Eg v = 0.4 mag and 0.65 
mag and varying widely (1 K kms -1 < Wco< 6-12 K kms -1 ) 
with Eg v in one of two narrow ranges, see Fig. 3: the same large 
scatter in CO column density with reddening and N(H) which is 
seen globally in galactic absorption line surveys also occurs in 
the single region studied here. The lower-reddening branch of 
the emission is spatially segregated to the galactic northeast of 
the star and has somewhat smaller mean integrated brightness 
Wco (Fig. 4) and velocity dispersion <x v (0.35 km s -1 vs. 0.42 
kms , Fig. 6). However, peak temperatures 8-12 K are present 
in both branches. 

The two most striking observational results of this study are 
the strong lines which emanate from the diffuse gas, up to Wco 
= 12 K km s~' , and the velocity structure present in the strongly 
emitting 12 CO line cores there (Figs. 7 and 8). The turbulent 



flows in this gas are in general directly visible as the spatially 
and kinematically resolved velocity gradients in simple, narrow, 
bright lines whose widths (cr v < 0.6 kms -1 ) are subsonic at dif- 
fuse cloud temperatures T K = 30 - 60 K (Fig 4 and 8). Other, 
locally-broader line profiles will likely be resolved into such ve- 
locity shifts of narrow line cores with higher (than 4') resolution 
although this remains to be tested. Conversely, it also remains 
to be seen whether CO line profiles which are subsonic but still 
supra-thermal at 2.7' resolution will show spatial or velocity sub- 
structure when mapped at higher resolution. 

We briefly discussed some modelling results of the formation 
and excitation of CO in diffuse media (Fig. 9). The brightness 
of the strongest CO lines can be understood by noting that the 
H2 column densities in the gas around f Oph are near the point 
(NCH2) ~ 10 21 crrT 2 ) where carbon fully recombines to CO at 
even modest densities n(H) = 200-500 crrT 3 , so that substantial 
columns of CO may form in gas which is at typical diffuse cloud 
temperatures (above 30 K). In turn, such densities will excite 
CO to the required degree even though they are far too low to 
thermalize the lower rotational level populations. 

In the range Wco ^ 10 K kms -1 it is the case that Wco x 
N(CO), as a consequence of very sub-thermal excitation. The 
CO-H2 conversion factor therefore varies widely in diffuse gas 
(because N(CO) varies rapidly with NCH2) and with great scat- 
ter), but it takes on values N(H 2 )/W C o = 2 x 10 20 H 2 cm" 2 /(K 
kms -1 ) in limited circumstances (Fig. 9 panel a at upper left) . 

If Wco x N(CO), a sky map of Wco like that in Fig. 1 is 
a map of the interstellar diffuse cloud chemistry. This should be 
contrasted with the more usual assumption of a constant CO con- 
version factor N(H2)/Wco, in which case a CO map traces the 
contours of the mass (i.e. the bulk molecular material) largely in- 
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dependent of physical conditions. Even the darker gas seen in the 
L204 complex South of f Oph is not immune to the influences 
of chemistry, which are clearly visible in the displacements be- 
tween regions of higher Wco and Eg_ v . 

The wealth of structure seen in the foreground CO bright- 
ness map has important implications for absorption line study 
of diffuse clouds. Since the viewing geometry is an accident of 
the relative locations of the Sun and £ Oph, our line of sight 
to the star could equally well have intercepted any of the wide 
variety of profiles shown in Fig. 7. This belies our ability to in- 
fer the general properties of the host gas from studies along any 
single sightline, no matter how comprehensive; studies of indi- 
vidual absorption sightlines must be viewed demographically, 
as datapoints within a large family of possible outcomes, even in 
nominally similar conditions. Even beyond this, there are some 
obvious fundamental limits to the use of absorption lines to de- 
rive the properties of the intervening gas. In the present case only 
a map could correct the false impression that the star is occulted 
by two separate foreground clouds. Likewise, the turbulent flows 
in the foreground gas appear clearly in maps of the gas but not 
at most individual positions, thereby conveying the false impres- 
sion of an overly-quiescent medium. 

We intend to map with higher angular resolution some re- 
gions of the L121 gas whose line profiles are thermal at the 2.7' 
NANTEN resolution, to see what kind of substructure might be 
present when some forms of line-broadening are absent. We will 
also map some L121 gas whose profiles are broader and whose 
velocity gradients are not fully resolved in Fig. 7 at 4' (0.16 pc) 
beam-spacing. 



This is the second paper in a loose series (see Pety et al. 
(2008)) which will also report observations of similar kinemat- 
ics on smaller angular scales and at higher angular resolution 
6"- 22" in other diffuse clouds of unknown distance whose pres- 
ence was first manifested in our mm-wave absorption studies of 
polyatomic molecules. In a subsequent paper we will discuss the 
geometrical and physical interpretation of the internal structures 
responsible for the flows seen there and in Fig. 7. 



— . 1997, A&A, 322, 962 

— . 2007a, A&A, 476, 291 

— . 2007b, A&A, 461, 205 

Liszt, H. S. & Lucas, R. 1998, A&A, 339, 561 

Maier, J. P., Lakin, N. M., Walker, G. A. H., & Bohlender, D. A. 2001, ApJ, 553, 
267 

Morton, D. C. 1975, ApJ, 197, 85 

Pety, J. & Falgarone, E. 2003, A&A, 412, 417 

Pety, J., Lucas, R., & Liszt, H. 2008, A&A, in press 

Rachford, B. L., Snow, T. P., Tumlinson, J., Shull, J. M., Blair, W. P., Ferlet, 
R., Friedman, S. D., Gry, C, Jenkins, E. B., Morton, D. C, Savage, B. D., 
Sonnentrucker, P., Vidal-Madjar, A., Welty, D. E., & York, D. G. 2002, ApJ, 
577, 221 

Savage, B. D., Drake, J. R, Budich, W., & Bohlin, R. C. 1977, ApJ, 216, 291 
Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Sonnentrucker, P., Welty, D. E., Thorburn, J. A., & York, D. G. 2007, Astrophys. 

J., Suppl. Ser., 168, 58 
Tachihara, K., Abe, R., Onishi, T., Mizuno, A., & Fukui, Y. 2000, Publ. Astron. 

Soc. Jpn., 52, 1147 
Van Dishoeck, E. F. & Black, J. H. 1986, Astrophys. J., Suppl. Ser., 62, 109 
— . 1988, ApJ, 334, 771 

Wannier, P., Penprase, B. E., & Andersson, B.-G. 1997, ApJ, 487, L165 

Wannier, P. G., Penzias, A. A., & Jenkins, E. B. 1982, ApJ, 254, 100 

Wolfire, M. G., Hollenbach, D., McKee, C. E, Tielens, A. G. G. M., & Bakes, E. 

L. O. 1995, ApJ, 443, 152 
Wolfire, M. G., McKee, C. E, Hollenbach, D., & Tielens, A. G. G. M. 2003, ApJ, 

587, 278 

Wood, K., Haffner, L. M., Reynolds, R. J., Mathis, J. S., & Madsen, G. 2005, 
ApJ, 633, 295 



Acknowledgements. The National Radio Astronomy Observatory is operated 
by Associated Universites, Inc. under a cooperative agreement with the US 
National Science Foundation. IRAM is operated by CNRS (France), the MPG 
(Germany) and the IGN (Spain). This research made use of the Simbad astro- 
nomical database and the NASA ADS astrophysical database system. This work 
profited from discussions of CO excitation with Michel Guelin. 



References 

Black, J. H. & Dalgarno, A. 1977, Astrophys. J., Suppl. Ser, 34, 405 
Burgh, E. B., France, K., & McCandliss, S. R. 2007, ApJ, 658, 446 
Cardelli, J. A., Mathis, J. S., Ebbets, D. C., & Savage, B. D. 1993, ApJ, 402, L17 
Crutcher, R. M. 1979, ApJ, 231, L151 

Dobashi, K., Uehara, H., Kandori, R., Sakurai, T., Kaiden, M., Umemoto, T., & 

Sato, F. 2005, Publ. Astron. Soc. Jpn., 57, SI 
Falgarone, E., Panis, J.-E, Heithausen, A., Perault, M., Stutzki, J., Puget, J.-L., 

& Bensch, F. 1998, A&A, 33 1 , 669 
Falgarone, E. & Phillips, T. G. 1996, ApJ, 472, 191 

Gaustad, J. E., McCullough, P. R., Rosing, W.. & Van Buren, D. 2001, Publ. 

Astron. Soc. Pac, 113, 1326 
Goldreich, P. & Kwan, J. 1974, ApJ, 189, 441 

Hartmann, D. & Burton, W. B. 1997, Atlas of galactic neutral hydrogen 

(Cambridge; New York: Cambridge University Press) 
Herbig, G. H. 1968, Zeitschrift fur Astrophysik, 68, 243 
Knapp, G. R. & Jura, M. 1976, ApJ, 209, 782 

Kopp, M., Gerin, M., Roueff, E., & Le Bourlot, J. 1996, A&A, 305, 558 
Lambert, D. L., Sheffer, Y, Gilliland, R. L., & Federman, S. R. 1994, ApJ, 420, 
756 

Langer, W. D., Glassgold, A. E., & Wilson, R. W. 1987, ApJ, 322, 450 
Liszt, H., Lucas, R., & Pety, J. 2006, A&A, 448, 253 
Liszt, H. S. 1979, ApJ, 233, L147 



Liszt, Pety and Tachihara: CO, reddening and turbulence around ( Oph, Online Material p 1 



Online Material 



Liszt, Pety and Tachihara: CO, reddening and turbulence around ( Oph, Online Material p 2 



1 1 1 1 1 1 1 1 1 1 1 in Hfi'.'.r.n-h 




0.2 



-i i 1 1 1 1 1 i.i 1 1 1 1 1 1 1 u 1 1 i ir 



0.0 



°- 5 « 1j0 

_ 00 



1.5 



[mag] 



2.0 



et al 



Fig. 10. Comparison of the limiting extinction A" of Dobashi 
( 2005| l with the limiting reddening data Eg_ v ofKcHegeT 
(T99S|. There is a left-offset of the A~ data of /8 mag, 
which is also seen toward ( Oph itself, see Sect. 6 of the text 
(available online). 
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7. Comparison of measured Eg v and A~ 

As an alternative source of extinction data we considered the re- 
sults of Dobashi et al. (2005 ) who constructed a sky map of limit- 
ing extinction A™ at latitudes \b\ < 40° based on star counts. We 
generated A" from a FITS file downloaded from the survey web- 
site, with pixels on a 2' grid. The features in a map of that data 
strongly resemble those shown in Fig. 1 but Ay = 0.25 mag to- 
ward £ Oph, which substantially underestimates the known fore- 
ground extinction since it is accepted that Eb-v = 0.32 mag and 
A v ~ 3.1x0.32= 1.0 mag. 

A more general comparison with the limiting extinction of 



Schlegel et al. ( 1998 1 is shown in Fig. 10: we generated A" and 
E£_ v on a 4' gncTTor 3° < / < 9°, 18° < b < 25° as in Fig. 1. 
The Eg v -intercept (0.32 mag) and slope (0.41) indicate an off- 
set of 0.32 mag/0.41 = 0.78 mag in the A" dataset with respect 
to E™ v , which corresponds to the disparity in Eg v and A" val- 
ues toward the star itself or to the minimum Eg v = 0.23 mag 
over the region. We infer that a uniform foreground component 
of the extinction, which might have renormalized the star count, 
is absent in the Ay maps for this region. 
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Fig. 11. Maps of various quantities, available online. Top left and right, peak and integrated brightness. Middle: mean velocity and 
velocity dispersion. Bottom, the velocity gradient.. 



